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Black Holes generate a particular kind of environment dominated by an accretion flow which 
concentrates a magnetic field. The interplay of gravity and magnetism creates this paradoxical 
situation where relativistic ejection is allowed and consequently high energy phenomena take place. 
Therefore Black Holes, which are very likely at the origin of powerfull astrophysical phenomena 
such as AGNs, micro- quasars and GRBs where relativistic ejections are observed, are at the heart 
of high energy astrophysics. The combination of General Relativity and Magneto-HydroDynamics 
(MHD) makes theory difficult; however great pionneers opened beautiful tracks in the seventies and 
left important problems to be solved for the next decades. These lectures will present the status 
, of these issues. They have a tutorial aspect together with a critical review aspect and contain also 

■ some new issues. Most of these lectures has been presented at the school on "Black Holes in the 

(~) Universe" at Cargese, in May 2003. 
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I. INTRODUCTION 



The existence of Black Holes is now very likely and established for some interesting cases: at the Galactic Centre, 
where not only the mass is known but also the spin; in the nucleus of some Seyfert galaxies where either maser 
effects have been observed in the accretion disk or where iron lines have been analysed and display a strong Doppler 
effect, that reveals a large mass concentration. The physics of Black Hole environments is getting more and more 
grounded and receives constant progress thanks to new generation of instruments, especially in X-ray and Gamma-ray 
' astronomy. 

The generation of relativistic jets is probably one of the most fascinating phenomenon attributed to Black Hole 
environments, where the accretion process governs most of the phenomenology. The fact that the most energetic 
phenomena in the Universe come from plasma ejections in the vicinity of Black Holes is a strong paradox that 
stimulates challenging theoretical investigations. 

Moreover, these relativistic ejections generate high energy phenomena, that are at the forefront of the astroparticles 
physics, which motivates the development of a new generation of instruments that flourish everywhere in the World. 

These lectures will introduce all the physics useful to built an understanding of what is happening with the rel- 
ativistic jets in Active Galactic Nuclei (hereafter AGNs) and micro-quasar jets, as well as in Gamma-Ray Bursts 
(hereafter GRBs). These topics incite the development of the physics of relativistic plasmas and relativistic MHD, 
which are one of the major front of modern plasma physics. 

These astrophysics topics are recurre ntly reviewed (| @ 0). However all the useful background of knowledge is 
briefly recalled, in order that the presentation be self-contained and read by a beginner. Furthermore, more detailed 
developments will be oriented towards numerous references, and some opened questions and trends will be presented. 
The content of the chapter is organised as follows. The main observational features of the relativistic ejections (AGNs, 
microQuasars and GRBs) are reviewed in section (2), with a presentation of the order of magnitudes of the main 
parameters, and of their performances as high energy engines. The physical properties of spinning Black Holes and 
the electromagnetism in Kerr metric is briefly recalled in section (3) in order to prepare the reader to the physics 
of jet generation. More detailed presentations can be found in text books (see for instance [5(j, HH, [zH)- Section 
(4) is devoted to the theoretical description of jets, including the presentation of the Blandford-Znajek mechanism of 
jet generation from a Kerr Black Hole. The properties of relativistic jets are different from those of non-relativistic 
jets, and they are emphasised in section (5), notably the difference in their confinement properties, the issue of the 
Compton drag suffered by relativistic jets, and the issue of the baryonic pollution of a pair jet. The other mechanism 
of jet generation from an accretion disk is presented in section (6), with an emphasis on the conditions, especially 
on the transport coefficients, to get launching. A presentation of the "two-flow" model is given in section (7) which 
consists in combining both a mildly or subrelativistic jet with relativistic ejections in the core, the subrelativistic 
jet playing a crucial role to circumvent the main difficulties encountered in the theoretical attempt to describe the 
generation of relativistic jets. The status of this field of research and of its new trends are summarised as a conclusion. 
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II. RELATIVISTIC EJECTIONS 



A. AGN Jets 



Some Active Galactic Nuclei (AGNs) exhibits powerful collimated jets. Although very menority, they are very 
interesting because, first the ejection phenomenon from a Black Hole environment is very intriguing, and also because 
jets are sources of high energy processes. There are two classes of jets, namely FR1 and FR2 jets (the classification is 
due to Fanarov -Riley (see Q))- FR2 jets are the most powerful, very collimated and ending with hot spots, whereas 
FR1 are weaker, shorter, displaying wiggles and no hot spot. 

From radio observations and Chandra measurements, the parameters of some FR2 jet hot spots are rather well 
estimated, such as Cygnus A and Pictor A hot spots (Wilson et al. 88]). Consider Cygnus A; the size (Rhs — 3 kpc) 
and the advance speed (Vhs — 0.1c) are known for a long time, the radio synchrotron emission has been mapped 
by VLA with accurate isophotes and the synchrotron self-Compton emission (SSC) is now measured in X-rays by 
Chandra. Thus the magnetic field intensity B can be estimated by the ratio of the luminosities, Lgsc 

over Lsyn ] and 

B ~ 2 x 10~ 4 G, a value close to the equipartition estimate. Assuming that the hot spots are powered by a mildly 
relativistic jet (/3jC), the jet power can be estimated and 

Pj ~ 6/3j x I0 45 erg/s . (1) 

This number is very comparable with the accretion power around a Black Hole of a few 10 8 A/ Q and thus indicates 
that a sizable fraction of the accretion power actually goes into the jets generation. 

In fact, as will be analysed later on in section (6), when an accretion disk launches jets, a sizeable fraction of the 
accretion power goes in the jets, and the bulk velocity essentially depends on the fraction of mass flux that is rerouted 
from the accretion flow to the outflow. 

Some population of AGNs exhibits high energy phenomena, such as hard X-ray emission and Gamma-ray emission, 
sometimes up to TeV range, together with a fast variability. These emissions come from a relativistic jet pointing 
towards the observer, as will be detailed in the subsection (2.3). Those sources are called "Blazars", which includes 
two categories of objects: flat spectrum radio quasars (FSRQ) and BL-Lacs, themselves divided in high BL-Lacs and 
low BL-Lacs. 



B. Eddington scaling 

Given the mass M of the black hole, several units of this physics can be defined (@). The gravitational radius 
r g = GM/c 2 is the basic length unit, from which an obvious time scale stems, namely r g = r g /c. The last stable 
orbit is at 6r g for a Schwarzschild Black Hole, at r g for a Kerr Black Hole at maximum rotation. The environment is 
made of a hot plasma undergoing an accretion flow and sometimes an ejection flow, both manisfesting through their 
radiation. The powers associated with these flows and radiation tightly depend on the accretion rate M a . Through the 
accretion flow, matter progressively goes deeper and deeper in an effective potential well that describes the effects of 
both gravitation and rotation. In Newtonian gravity, this is simply $ e /y(r) = —GM/r + l 2 /2r 2 where I is the specific 
angular momentum of the matter element. The effective potential is deeply modified by General Relativity effects, 
described by the two Black Hole parameters, its mass M and its reduced spin parameter a s (\a s \ < 1). Accretion 
power characterises the rate of energy liberated by matter flowing from infinity down to the last stable circular orbit; 
below the last stable orbit, matter falls almost freely towards the Black Hole horizon, without significantly radiating. 
In Newtonian dynamics, the accretion power is 

P. = \M„™ . (2, 
2 r lso 

More generally, it is given by 

P a = e a M a c 2 , (3) 

where e a is called the accretion efficiency and includes relativiqstic effects; its value is 6% for a Schwarzschild Black 
Hole, 42% for a Kerr Black at maximum rotation. 

The Eddington luminosity is a useful reference for the accretion power: 

L<L E = —GMm p ~ 1.3 x 10 46 (— ^— )erg S /s (4) 
<7t 1U°A2© 



3 



For a given efficiency, it defines also a reference for the accretion rate: M a = mMEd an d P a — ttiLes,- 

The accretion power is shared between the radiation flux of the disk and the power in the jets it launches. L = xP a 
and Pj = |(1 -x) -Pa- 

The luminosity emitted by the disk depends on its opacity, the most powerful emission occurs when the disk is 
optically thick. In the hottest part of the disk and of the jet, the Compton opacity is the relevant one: r = cr-rroro. 
In standard accretion disk (hereafter SAD) r > 1 ([fH), whereas weakly radiative tori with r < 1 has also been 
considered, as, for instance, advective dominated accretion flow (hereafter ADAF) (|H^). 

The black body temperature of the SAD is maximum at a few last stable orbit radius, 8r g typically (for a s = 0) 
and determines the typical frequency range of thermal emission. 



For the Eddington luminosity, 



M 

10 ^wsj*"- < 6 > 



The corresponding maximum pressure is P* = ^a s jT^ ~ 3 x 10~ 2 L/(7rr 2 c). 

The magnetic field plays a key role together with gravitation in the physics of accretion, ejection and its associated 
high energy manifestation. It is concentrated by the accretion flow; Dynamo action could also be at work in accretion 
disk f |l5j1. For physical reasons, it is convenient to compare its intensity with the value obtained by assuming 
cquipartition of its pressure with the kinetic pressure (of the plasma or the photons) on the equatorial plane of the 
accretion flow (an excess of magnetic field close to the central black hole would trigger interchange instabilities that 
would expel the excess). For a Schwarzschild Black Hole, the equipartition magnetic field for Eddington luminosity 
is thus 

B -^-^^Jk>-'»G. (7) 

In the case of optically thin tori (which is likely the accretion regime close to a spinning black hole Q [24|). the 
magnetic field pressure does not exceed a fraction of pf2 2 r 2 , and, because the mass density p is much weaker than 
in the case of an optically thick disk (naxh -C 1 with h ~ r), the magnetic field intensity can be stronger than the 
previous estimate, but not more than a factor ten, and follows the same scaling with M (however its dependence 
with a s need some further examination). The electric current flowing in the accretion-ejection system is particular 
important as will be seen further on; (B*r g ) can be considered as its fiducial value. The magnetic field and the current 
together with rotation allow to develop a power for the generation of jets. The magnetic field plays an important role 
in the transport phenomena (diffusion, effective viscosity, thermal conduction etc.) and in the processes of particle 
acceleration. With simple arguments, a fiducial number can be proposed for the maximum energy of the Cosmic 
Rays accelerated in the jets. MHD perturbations (waves, reconnections and shocks) in the jets can no more accelerate 
particles whose Larmor radii would exceed the radius of the jets; the following maximum energy of Cosmic Rays is 
thus obtained: e ma x — 10 20 TZ(BR) eV (with B measured in Gauss and R in pc); and BRj ~ ( 10 sm ) 1 ^ 2 G x pc 

e max ~ io 20 rz(-^-) 1 / 2 e y 

V 10 8 Mq ' 

As a result of particle acceleration, all kinds of high energy radiations could be emitted by the jets: synchrotron 
emission by relativistic electrons from radio band up to IR, UV, and sometimes X-rays ranges, inverse Compton 
emission in gamma ray range, sometimes up to TeV range, p7-Neutrinos emission as a signature of very high energy 
cosmic rays. Black Hole formation and jet acceleration generate gravitational waves that could hopefully be detected 
in the future f[68|). 

These sources of radiations are compact in a precise physical sense, and characterised by the "compactness param- 
eter" : 

W-7^7 = ^?- (8) 

The compactness is related to the opacity to 7-rays due to the process of pair creation by photon-photon interaction. 
The threshold of the process is such that €162 > 2/(1 — cos 9) where e = hv/m e c 2 . A radiation field of soft photons, of 
energies e s , with a spectrum S(e s ) oc ej a , constitutes an absorbing medium for gamma photon. For gamma-photons 
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of energies e, the opacity to pair creation is given by r 77 ~ 0.2£ c e a ( 73]). Thus the higher the energy of the 7-photon 
the stronger the absorption. The nuclei of AGN and of microquasars are clearly optically thick to gamma rays since 
the compactness is 10 2 — 10 3 and can be richly populated with e + e~-pairs. Compactness characterises not only the 
opacity to gamma-rays, but also measures the strength of the Compton drag exerted on a relativistic flow, as will be 
seen later on in subsection (5.1) and section (7). 



C. Gamma- Ray Emission and Relativistic Kinematics 



Gamma rays can escape only at remote distance and thanks to relativistic motions. A source of photons in 
relativistic motions (/3o close to 1) with a bulk Lorentz factor T = (1 — /3o) -1 ^ 2 > 2, is a cause of a Doppler shift, of 
a Doppler beaming, of an apparent superluminal velocity and a shorter apparent variability. 

• The photons energy are blue shifted when the flow points toward the observer at an angle 6: 

v' = T(1- 0co&9)v = v/6 . (9) 
Synchrotron and Compton peaks are blue shifted v \, s — Sv com . 

• The variability time scale appears much shorter than in the co-moving flow: At b s — At com /S. 

• The emission is beamed. Because I{v, n)/V 3 is proportional to the number of photons in the state defined by v 
and n, it is a Lorentz invariant. Thus 

I{v,n) =5 3 I com (p/S) . (10) 



• The velocity of the source can appear to be faster than the speed of light. Indeed, during St, the source 
travels a distance vSt sin (9 as projected on the sky plane; whereas the apparent duration of the travel is St app = 
(1 — fa cos 9)St. The apparent velocity is thus 

sin 9 

Vapp = v- , 11 

1 — po cos 9 

and reaches a maximum Tv for cos 9 = [3q or sin 9 = 1/T. 

Moreover, if the source of observed photons is moving relativistically, the estimate of the intrinsic compactness 
must be rescaled compared to the isotropic assumption made for calculating the luminosity (L' — S^ i L, E' = SE, 
At' = SAt bs) so that £ c — 6~ 5 £ c \ b s - For a = 1, t 77 = <5~ 6 r 77 | fc s . For AGN, observation of GeV photons requires 
S > 5. For GRBs, the compactness is very high (£ c ~ 10 12 ) and observation of lOMeV photons requires S > 100. 
The opacity to pair creation leads to an observed break in gamma-spectra in MeV range, at 

e b = S(0.2 4r 1/Q • 

These relativity kinematics effects are crucial to analyse the Blazar spectra with their double peaks, synchrotron peak 
in X-rays and SSC peak in gamma-ray, together with variability (see |38l ] for recent fits). The observation of gamma 
rays is a strong argument in favour of relativistic jets. This argument adds to the historical argument inferred by 
Martin Rees in 1966 (|z3) m favour of relativistic motions in AGNs in order to prevent the Compton catastrophe 
(0)- 



D. Microquasars 

Galactic sources, which have been discovered during the last decade 0] H^, have features similar to radio quasars 
: a bright nucleus with twin radio jets displaying superluminal motions and hard X-ray IC-emission; they have been 
called micro-quasars. 

Most of the phenomenology can be understood by simply rescaling the Eddington fiducial numbers previously 
introduced that depends on the Black Hole mass only: 10 7 — 10 9 M Q for AGNs and ~ IOMq for micro-quasars. 
The inverse Compton process on accretion disk photons sensitively changes in the case of micro-quasars because 
the disk black body radiation is in the soft X-rays and thus the inverse Compton process is deeply modified by the 
Klein-Nishina regime (see subsection (7.2)). 
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At first look, the proximity of micro-quasars seems interesting, letting think that they will be studied more easily 
than remote AGNs. In fact, in terms of gravitational radii, the angular resolution is much better with AGNs at a few 
Gpc than with the closer micro-quasars. The advantage of the smallness and closeness of micro-quasars is the much 
shorter time evolution. Indeed, connection between the accretion process and the ejection has been nicely revealed by 
the changes of regime (|23|)(|7(). These variations are very long compared to the Keplerian time at the inner orbits of 
the accretion disk, and similar variations could not be observed in AGNs during a human life. On the other hand, the 
shortest variability at Keplerian time is better recordered in AGNs, especially in Gamma-Rays; next X-ray satellites 
(as for instance SIMBOL-X) will probably be able to follow millisecond variations in micro-quasars. 

As source of cosmic rays of energies up to 10 16 eV, a neutrinos emission resulting from jr^-collisions could be detected 



E. Gamma Ray Bursts 



Two kinds of progenitors have been considered to give rise to the GRB phenomena: coalescence of compact objects 
and collapse of supermassive stars. Besides two kinds of origination, there are two classes of models for the description 
of the GRB phenomenon itself : the "fireball" model (|z3]) and the electromagnetic model (0). 

The fireball model is based on the assumption of an initial state of high relativistic enthalpy that gives rise to 
a relativistic flow through adiabatic expansion. The energy involved in the explosion is E ~ 10 51 erg; the initial 
temperature is therefore To ~ 5MeV. The dynamics depends on the baryon load of the relativistic wind; it is 
measured by the parameter r/ = E/M^c . The bulk Lorentz factor saturates at the value r\ after a travel from few 
gravitational radii (ro) to r s — rjro. Then the ejected shells start broadening at r& = rj tq. Nowadays, the value 
7] ~ 300 allows the best fit of most of the observed GRBs, but no physical argument to select this value has been 
proposed so far. 

To explain the fast irregularities of the light curves of many GRBs, the fire ball model has been remodeled by 
considering multiple shells instead of a single one [44|. The date of collision between two shells launched at a time 
interval At is given by 



r 2 r 2 

p2 _ p2' 
1 2 1 1 



*c = 2t^2 At • (I 2 ) 



The maximum duration of the prompt emission is At max ~ 2T 2 At w , where At w is the duration of the relativistic 
ejection. The observed duration of the prompt emission is At f, s = (1 — /3)At max ~ At max /2T 2 ~ At w ; therefore the 
observations give access to the duration of the ejection, and two categories of GRBs have been observed: short ones 
with < At w >~ 0.1s and long ones with < At w >~ 10 s. Similarly, the observed variability time corresponds to 
irregularities at the origin of the flow with a similar time scale, At b s ~ At var . Millisecond variations are sometimes 
observed; they correspond to shell collisions at a few r\, = r/ 2 ro. The proper length of the flow is Iq — f3cTAt w which 
indicates that the jet length is smaller than the distance where the shells start to collide each other {j3cAt w < r^). 
The flow and the magnetic field are likely disconnected from the central engine when the prompt emission starts. 
Fairly collimated (for a few GRBs, the average solid angle is estimated as Q,/Ait ~ 2 x 10 -3 ), these relativistic flow 
can be considered as transient jets. However there is no physical argument either for explaining the duration of the 
flow. They end in the form of a strong relativistic shock with the interstellar medium, or the ejecta medium, that 
produces the so-called "afterglow", which has been well studied and brought support to the model. 

The equipartition magnetic field can be as high as B ~ 10 15 G at the initial state. The magnetic pressure decreases 
in the expansion like the kinetic pressure, namely P m cx P oc r~ 4 , as long as the magnetic field is considered to be 
tangled, carried by the shells and disconnected from the source (j45j). In the canonical fireball model, the magnetic 
field can never dominate. 

In the electromagnetic model, a "magnetar" object is initially created ([4lj). whose environment contains a rather 
cool plasma highly dominated by the magnetic field, and which is in fast rotation. Therefore the relativistic wind is 
no more due to adiabatic expansion, but to a powerful Poynting flux. This model is presented in another chapter of 
this book (see Lyutikov's and Blandford's chapters). 

GRBs are popular events for high energy astrophysics because, not only they are very powerful gamma-ray emitters, 
but also they are the most promising sources of UHECRs,with an associated high energy Neutrinos emission (j87(), 
and probably detectable sources of Gravitational Waves ([63]). 
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F. Performances of the Ejections 

In the tabled the order of magnitude of the main parameters of the three kinds of relativistic sources are indicated: 
their power, the mass of the central black hole, the gravitational time that indicates the fastest variability time 
scale, the bulk Lorentz factor of the jet, and the equipartition value of the magnetic field in the vicinity of the black 
hole. Thanks to magnetic irregularities and fast motions, both processes of Fermi acceleration, by MHD turbulence 





Power (erg/s) 


M BH (Mq) 


Tg ( s ) 


r 


B (G) 


AGNs 


^q44-46 


w r-9 


1Q 2-4 


~ 10 


~ 10 4 


/^-quasars 


~ 10 38 


~ 10 


10~ 4 


~ 10 


~ 10 s 


GRBs 


10 49-52 


10-100 


~ 10~ 3 


- 300 


~ 10 15 



TABLE I: Performance of the ejections. 



and shocks, generate ultrarelativistic electrons that radiate synchrotron emission. Lorentz factors achieved by the 
electrons in the co-moving frame are on the order of 10 6 in AGN and micro-quasar jets. Therefore jets produce a 
synchrotron emission in Radio band and sometimes up to X-rays. The main gamma-ray emission of GRBs is produced 
by ultrarelativistic electrons that radiate synchrotron keV-photons in the co-moving frame. 

Acceleration of cosmic ray protons is expected in these three kinds of relativistic flows. As previously explained, 
their energy is necessarily smaller than the confinement limit which can reach a few 10 20 eU at the beginning of the 
AGN jets where the electric current is maximum. This concerns FR2 jets, for FR1 jets like Centaurus A do not have 
a sufficient product BRj (see Morganti et al. |5l|), thus the current is too low for achieving ultra high energies. The 
same argument for micro-quasars leads to 10 16 eV. AGN hot spots are also considered as possible sources of UHECRs; 
however this issue will not be discussed in the chapter. 

In GRBs, the limiting size for the maximum cosmic ray energy is the width of the shells at the beginning of the stage 
of internal shocks, namely at r&, where their width is AR ~ r/T in the co-moving frame. Thus e max ~ 10 20 Z(Br)b 
and since rt ~ 10 12 cm, the generation of UHECRs would require a magnetic field of order 10 6 G, and thus ~ 10 15 G 
at the initial state. It turns out that the terminal shock is not able to generate the expected flux of UHECRs, but 
this will not be discussed in this chapter either. 

III. SPINNING BLACK HOLES AND RELATIVISTIC MHD 

The main point of the jet formation is that opened magnetic field lines thread either an accretion disk (or torus) 
or a spinning black hole and the rotation of the field lines differs from the rotation of the matter or the black hole. 
Thus a poloidal current is driven whose associated toroidal field is responsible for 

• transfering the angular momentum from the rotator through the magnetic stress (and torque) —B^Bp/Aw; 

• generating of poloidal Poynting flux E x B^,/4tt, the electric field being generated by the rotation of the field 
lines; a poloidal force also pushes matter along the poloidal direction; cx — V(r£?0) 2 

• a confinement is possible because of the magnetic tension —B^/Air. 

The general pattern of the accretion-ejection flow is sketched on figure (1). The poloidal magnetic field (blue lines) 
has a bipolar configuration (a quadrupolar configuration is inefficient to extract the angular momentum). The poloidal 
current (yellow lines), associated with the toroidal field generated by differential rotation, flows towards the central 
object and the inner part of the accretion disk; the current closure is generally insured at remote distance from axis 
with current lines emerging from the disk and unfolding in the jets (see [2(J). Note that the electric field (red lines) 

FIG. 1: Sketch of the general pattern of the MHD flow. In blue the field lines, in yellow the current, in red the electric field, 
in green the matter flow. The dashed line suggests the "watershed" surface. 



points towards the axis almost everywhere except close to the Black Hole horizon as will explained later on. Part of 
the accretion flow (green lines) is rerouted by the magnetic field and supplies the jet, it crosses the magnetic surfaces 
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in the accretion dissipative region and then the outflow runs along the poloidal field lines in ideal MHD regime. Thus 
there is a "watershed" surface in the disk that separates the matter flow that will actually be accreted by the black 
hole and the matter outflow. Above the black hole, no matter would be expected; however pair creation takes place 
there and the watershed surface separates the pair flow that falls onto the black hole and the pair flow that is expelled 
along the axis. 

A. Summary of the "3+1" split of the Kerr metric 

The Kerr metrics of space-time is the solution of General Relativisty equations (see with two symmetries: 
stationarity and axial symmetry. Whereas the Schwarzschild solution has a point singularity at r = 0, the Kerr 
solution has a ring singularity at r — and 9 = ir/2. It depends on two parameters only: the mass M (or r g = GM/c 2 ) 
and the spin parameter a s . The angular momentum of the Black Hole is J = a s r g Mc. In astrophysics, the third 
parameter of black hole solutions, namely the electric charge, is always zero, in practice. The Boyer-Lindquist system 
of coordinates for a spinning black hole is a suitable universal frame for the space-time 4- manifold £4. 

ds 2 = -a 2 dt 2 + g v} (dx l + (3 i dt)(dx j + (3 1 dt) , (13) 

where dt stands for cdt. The main relativistic effects are described by the two following functions of r and 9: the 
"lapse" function a = h 1 h 2 /h 3 (multiplied by c 2 in usual units) and the GM-potential w — {2a s /h 2 ){r/r g ) (multiplied 
by c/r g in usual units). The "frame drag", specific to Black Hole rotation, is characterised by (3$ = —u (j3 r — (3 e = 0). 
The metric elements are 

9rr = (h 2 /hi) 2 , g gg = r 2 h\ , g^ = r 2 = r 2 (h 3 /h 2 ) 2 sin 2 9 , 
where the three functions hi, h 2 and ft. 3 depend on r/r g , and a s . 

r 9 r 9 

h 2 2 = (-) 2 +a 2 s cos 2 e 
r g 

hi = C-^ + a 2 ) 2 -a 2 hi sin 2 6 
r g 

The BL-system of coordinates displays a universal time t that allows to slice the four dimensional space-time into 3D- 
subspaces at constant t, these subspaces being isomorph to the Riemanian 3-manifold T 3 with metric dl 2 = gtjdx 1 dx^ , 
endowed with spatial curvature that characterises gravitation and that prevents the extension of any local inertial 
system of coordinates to a global one. At any point of a "Fiducial Observer" (FIDO) is at rest with respect of 
the local time r and thus moves with respect to the BL-system; it rotates at the sheared angular velocity oj. The 
conversion of the local time r such that ds 2 = dr 2 — dl 2 to the universal one is characterised by the lapse function 
a = dr/dt that measures the gravitational red-shift. The temporal curvature makes the global extension of the local 
time impossible. Note that gtt = —a 2 + w 2 f 2 = — {h\ — a 2 sin 2 0)/h 2 and gt<j> — g<f,t — —lot 2 . 

The length of the circumference defined by the set of points of given (r, 9) is 2irf ; and f is called the " circumferential" 
radius and reduces to rsin# at remote distance fro m the B lack Hole. The horizon of the black hole is such that 
a(r, 9) = 0; this is a 2-sphere of radius rn — r g (l + \J\ — a 2 ). The area of the horizon is A — ^{r 2 H + a 2 r 2 ). 

If one looks for FIDO that would be static, namely, at rest with respect to the BL-system, then, because of 
Black H ole rotation, it turns out that none exists inside the so-called " ergosphere" , a 2-manifold such that r = 
r g (l + \J\ — a 2 cos 2 9) (flattened by rotation effect), which corresponds to g u = 0. They are constrained to rotate 
with an angular velocity between f2i and where £l\, 2 = lu ± a/f. It is remarkable that these two angular velocities 
merge at the horizon: Qi = £^2 = w \h- That angular velocity u) \h is uniform over the horizon ("solid rotation") and 
can be defined as the Black Hole angular velocity £Ih', its relation to the spin parameter is given by Q h — a s c/2rH- 
FIDOs inside the ergosphere necessarly rotate and for the angular velocity to they have no angular momentum with 
respect to the BL-system; so they are also called ZAMOs (Zero Angular Momentum Observers). 

Because of stationarity with respect to the universal time, an energy invariant e exists and because of axisymmetry, 
an angular momentum invariant L z exists for geodesic motions. In terms of Killing generators of the symmetries V M , 
the invariants are p^V^ 1 . With respect to the BL-reference frame, 

e = -po = -(gttp a + gt^p*) 
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and 

L z = p,/, = g^ + g^tp 
Thus, combining both relations, the energy can be written: 

e = ~Po = a 2 p° + u)L z 

It is interesting to write these relations in terms of FIDOs measurement of energy and angular momentum. The 
generator of axisymmetry is unchanged: d<p; consequently the local measurement of the angular momentum gives 
the correct invariant L z = p^\fido- This is different with energy that varies locally; and the generator of temporal 
invariance <9 t = ad T — ojd<f,. Therefore — po = — ccPqIfido + ^>P4>\fido- Thus the energy invariant, corresponding to 
the energy measured at infinity, is given by 

£oo = cte + ujL z . 

In the ergosphere, it is possible to get a negative energy for L z < 0. This is at the origin of the Penrose process. Let 
a body in a state of positive energy eo in the ergosphere that decays into one piece with positive energy t\ following 
an outward geodesic and another piece falling down with a state of negative energy ti- Then e± = eg — 62 > eo, thus 
the outgoing body emerges with more energy to the expense of the rotation energy of the Black Hole. 

In the "3+l"-split formulation, the forces exerted on matter are measured by FIDOs in Tz (see [7^). The gravita- 
tional interaction is described by a local gravity field g = — V log a (times c 2 in usual units) and by a " gravito-magnetic" 
(GM) field Hij, defined by ffy = — Vj/3j, which acts on a particle like a magnetic field: dpi/dr = H^p 3 . Indeed, like 
the magnetic field that complements the electric field, as required by relativity, when charged particles are moving, 
similarly the GM-field complements the gravity field when massive particles are moving. The frame angular velocity 10 
is therefore called "gravito-magnetic" potential. The GM-field is responsible for a precession effect, the Lense-Thirring 
effect. Moreover, it acts on a viscous disk in such a way that matter is enforced to rotates in the equatorial plane 
orthogonal to the Black Hole rotation axis. Spinning black holes maintain a stable direction for jets that are supposed 
to be launched along the rotation axis; and if the alignement of the disk axis and the Black Hole axis is not completely 
achieved, the Lense-Thirring effect would produce an observable precession of the jets. 



B. The Membrane Paradigm 

An important issue about FIDOs is that, for BL-frame, they run faster and faster against the geodesic flow when 
closer to the horizon and all the motions seem to synchronise with FIDOs motions, because the velocity is such that 
\v\ — ct" 1 !^ + 0\ < lj which implies |4| + /3\ as small as a. Deviations to FIDOs motions close to the horizon are 
very slow motions. 

Transverse electric and magnetic fields diverge like a -1 close to the horizon and the most diverging contributions 
satisfy the property of a converging transverse electromagnetic wave in vacuum, namely B = —ri x E (|z3). In the 
membrane paradigm, boundary conditions are set up at a 2-manifold defined by a(r,0) — an <C 1. The normal 
electric field is accounted by an effective surface charge such that <jh = E n /Air and a vacuum surface resistivity 
(Rh = 4-7T = 377 Ohms) accounts for a Active current that closes the electric circuit such that the surface current Jh 
generates the magnetic field B = AttJh x n = —ri x E with E = RhJh- 

At the membrane it is convenient to rescale all the diverging fields by multiplying them by an f|79jp. In particular, 
one defines the local surface gravity gn — ci^g. 

An effective Ohmic dissipation is thus proposed by the membrane paradigm, which therefore implies a possible 
entropy production. Indeed hiding all the "sluggish" dynamics inside the membrane amounts to ignore a lot of 
informations and to constitute a state of high entropy The local quantum fluctuation time scale is such that St ~ c/ gn 
and a corresponding fluctuation temperature, Tjj ~ H/St, has been proposed by Hawking ('|33||'): 

rp rp h 9H ™ ~<" 2 

+ H = Oi H l = = 



2tt c 8ttM 



where mp is the Planck mass {rap = (fic/G) 1 / 2 ~ 10 19 GeF); this temperature is very low for astrophysical black 
holes and the quantum evaporation very inefficient. The entropy Sh is 1/4 times the ratio of the horizon surface area 
over the square of the Planck length (Bekenstein) [j| : 



S H = 2tt(1 + Vl-a 2 ) 



,.2 
9 
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which is an extremely large number for astrophysical black holes [lp = (foG/c 3 ) 1 / 2 ~ 1.6 x 10~ 33 cm). Although 
ineffective in astrophysical applications, these concepts are very interesting and profound, and they are very relevant 
in string theory at Planck scale. This entropy has recently been derived from quantum gravity theory (|z3) by 
counting the number of quantum states compatible with the macroscopic parameters of a black hole. 



C. Electromagnetism in Kerr metric 



In FIDO's space electromagnetism can be described by local Maxwell equations with local time. Then the 
universal time can be introduced, which changes Maxwell equation by inserting the lapse function a and the GM- 
potential w at the appropriate places. The electromagnetic field is thus coupled with the gravitation field. The 
coupling with matter is insured by Ohm's law (i.e. E + v x B = rjj for non-relativistic flows). In most astrophysical 
applications, the plasma is very conducting and the resistivity effect (77) is negligible (the Reynolds magnetic number 
is very large); the magnetic field is thus frozen in the plasma flow and the electric field is just there to compensate 
the effect of motions across the magnetic field, it vanishes in a frame comoving with the magnetic lines. 

In axially symmetric configuration, it is convenient for both analytical and numerical calculations to describe the 
electrodynamics with only two scalar functions a and b, that are defined as follows. Any vector field A can be splitted 
in two pieces, the poloidal one, A r e r + A e eg, and the toroidal one, A^e^. The basis vectors are not necessarly unitary; 
it is even convenient to use the natural basis associated with the coordinates and often labeled {d r ,dg,d^,} and the 
field components are the contravariant components. The covariant components are such that Aj — gjkA k . Covariant 
components change like components of a differential form under a change of coordinates, whereas contravariant 
components change like components of a first order differential operator; which justifies the notation with partial 
derivatives for the natural base of the vector fields; similarly the natural base of covariant field can suitably be written 
{dr, dO, d(f>}. The same when the time coordinate is introduced, one has a fourth base element dt for the contravariant 
field base and dt for the covariant field base. The poloidal magnetic field is derived from the covariant component 
of the toroidal vector potential; and one notes a = A^ — f 2 A^ — fA^ where A^ is the usual toroidal component for 
unitary basis vectors. The magnetic flux across an horizontal disk of constant r and 9 points is 2ira. 

B p = -^xVo (14) 



The poloidal component of Ohm's law reads: 



1 d 
act 



^-a + av j V i a = (16) 
at 



This clearly expresses the frozen in condition of the magnetic flux in the flow. Similarly the poloidal current is 
derived from the covariant component of the magnetic field B^, noted b. From Ampere theorem in magneto-static 
approximation, the current intensity crossing the same disk as before is / = 27ra6/47r = ab/2. The toroidal part of 
the Faraday induction equation reads: 

d 

—b + f 2 div[ab/f 2 v p - (aSl + w)B p ] = (17) 

The matter is supposed to rotate at an angular velocity SI = v^. 

Consider a stationary magnetosphere. Relativistic MHD will be reduced to five invariant quantities and a non- 
linear partial differential equation governing the transverse equilibrium of each magnetic surface, the so-called Grad- 
Shafranov equation (see [g). The magnetic surfaces are characterised by constant values of a. The toroidal electric 
field measured by FIDOs thus vanishes and the poloidal flow is along the poloidal field lines (v p x B p — 0). The 
matter flux, with a bulk Lorentz factor F can be written as follows apTv p = kB p and, since in stationary flow, it is 

divergence free like the magnetic field, B p .\7k = 0, and the function k, which is the flux of matter per magnetic flux, 
is constant along the field lines, it depends on a only; this is the first of the five invariants: 



(18) 
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Inserting the mass flux in term of the magnetic field in the induction equation l|17|) . it is easily found that the 
angular velocity 0„, defined by 

ft* = aQ + u>- -^6 , (19) 
pi r z 

is also constant over each magnetic surface (tt*(a)). It turns out that this angular velocity is the angular velocity 
of the frame in which the electric field vanishes, and one calls it the angular velocity of the field lines. Indeed the 
poloidal electric field is obtained from Ohm's law by inserting 

b _ n + \ 5 (20) 

K 

and one find 

-* fL — uj , 

E = Va . (21) 

ac 

A way to see that the field lines rotate despite the axisymmetry is to consider the toroidal Poynting flux which reads 

(£Ln^f£ll _e ■ it clearly shows that poloidal field energy rotates with respect to ZAMOs at an angular velocity (£!* — u)/a 
and O* in the BL-system of coordinates. Equation l|20|) expresses the current generation due to the shift of angular 
velocity across the field lines. Having got the electric field and the current, the main ingredients to estimate the 
Poynting flux are at hand. 

IV. GENERAL PROPERTIES OF JETS 

The general properties of jets depend on the transport of the angular momentum they carry from the central 
rotator, spinning Black Hole, accretion disk or torus. Their power, motions and collimation are tightly related to the 
transport of angular momentum which is itself linked to the electric current flowing in the jet. 

A. Transport of angular momentum 

The specific angular momentum is the toroidal covariant component of the specific momentum: I = u^c = Fv,/) = 
Tttf 2 . 

Its transport equation expresses the conservation of the global angular momentum of the system: 

— [(e + P)Tl + sJ ,,m - ) ] + div[(e + P)aTlv p - ^(bB p + E^E p )} = . (22) 

In stationary state (E$ = 0), inserting the mass flux in term of the poloidal magnetic field leads to a third constant 
£*(a) on each magnetic surface that expresses the conservation of the total angular momentum: 

ab 

wl-- — = Ua) , (23) 

47TK 

where w = (e + P)/pc 2 is the specific relativistic enthalpy. It is noteworthy to remark that in the region where the 
magnetic field dominates (quasi force free regime) the total specific angular momentum carried by matter is directly 
related to the current (J = ab/2) : Z*(a) ~ —I/2itk. For consistency the electric current must be a poloidal invariant 
in this regime, which is true for the force free approximation. Since 7 < because generated by positive rotation with 
respect to the direction of the poloidal field, the specific angular momentum has the sign of mass flux (k) . Infalling 
matter carries negative angular momentum towards the ergosphere of the Black Hole. 

B. The Alfven surface 



The Alfven surface is one of the most important critical surface of the flow. A fast flow is supposed to pass 
smoothly across this surface in order to prevent backward Alfven waves to accumulate there and generate a strong 
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current perturbation. Relativistic Alfven waves propagate at a velocity Va with corresponding Lorentz factor Fa such 
that 

4^ = ^- (24) 
e + P 

A local Alfven Mach number m can be defined such as m 2 = c2m / ; where involves the poloidal magnetic field 

only. The Alfven surface is such that m 2 (r,6) — 1. This criticality appears when the current and the angular 
momentum are derived from both induction equation (|2U[) and angular momentum equation l|23|) . Indeed the two 
following relations are deduced : 

-^-[m 2 - 1)6 = h - - (O* - uj)Tf 2 , (25) 
4ttk a 



w(m 2 - 1)1 = m 2 U - -(fi* - w)Ir 2 (26) 
a 

The regularity conditions imply that the magnetic surfaces must cross the Alfven surface at a distance from the axis 
such that 



^2 _ -2 _ 



r 



The values of 6 and I at the Alfven surface are then obtained as the ratio of the derivatives along the poloidal field of 
the vanishing factors (l'Hopital's rule). 

The important aspect of the Alfven surface is that it marks a transition between a regime of acceleration in widening 
magnetic surfaces and an asymptotic regime of quasi-ballistic motions, and also a transition in the transfer of angular 
momentum. Indeed, at the beginning of the flow where m 2 <C 1 (sub-Alfvenic regime) and f < fi, equation 12fi|) 
shows that I <C I* and the contribution of the magnetic torque to the angular momentum invariant Z* dominates. In 
the opposite super-Alfvenic case (m 2 3> 1), with f ^> ¥ai I — * I*, although the angular velocity decreases. In fact I is 
already close to Z» at the Alfven surface. 

A determination of the Alfven radius is crucial but it depends on technical details. A natural estimate can be found 
in non- relativistic regime, where it essentially depends on the mass flux Mj(a) inside a flux tube defined by a, on 
the field line rotation fi* and the magnetic flux 2ira. Assume Vl*fA of the order of the Alfven velocity at the Alfven 
surface, then, since the kinetic and magnetic energy density are equal at this surface, the energy fluxes are the same, 
and one can write: 

1 • B 2 TT 2 a 2 

-M (a)ttlr 2 A = \ Q -^n*r A Trr 2 A = A il, , (28) 

Z 47T r a 



where Aq ~ A ~ 1. Therefore the following estimate is derived 



r A = C-r^) V3 ■ ( 29 ) 



MjVL 



A relativistic estimate will be derived in a subsection 1IVDII 



C. Energy transport, Bernoulli invariant 

In stationary state, the conservation of energy is expressed by a divergence free flux: 

div(a 2 S) + ia 2 (V,/3 fe + V k ^)T^ k = (30) 

The multiplication of the local energy flux by a 2 is explained as follows: one factor comes from the time rescaling, 
another one comes from the conversion of the energy locally measured by FIDOs to the energy-at-infinity. The poloidal 
contribution of matter is 

§(m) = ( e + p )r 2^ . (31) 



12 



the contribution of the electromagnetic field is 

§{ e. m , = ExB, = (32) 

47T a 47T 

The second term of <|3(J[) is a non-trivial contribution of the gravito-magnetic field. It can be cast in the form of a 
divergence of the following field: 

a 2g(g.ra.) = Kl ^B p = ( WKl „ ^L) u B p . (33) 

Ait 

It is important to realise that, since b < and f2* < u> close to the Black Hole, the Poynting Flux is converging. But 
the contribution of the gravito-magnetic flux cancels the w-terms, so that the sum of electromagnetic and gravito- 
magnetic fluxes makes an outgoing energy flux. Which means that the gravito-magnetic effect is the one which makes 
an energy outflux from the rotating black hole. In particular, as long as the matter contribution is negligible (like in 
force free regime), the energy flux is 

a 2g(e.g.m.) = _^_^B p = ^-^B p . (34) 
47T ZTT 

However as one looks at a remote distance from the black hole, the GM-potential vanishes and the outgoing energy 
flux becomes essentially the Poynting flux. The form of the GM-energy flux, involving the flux of angular momentum 
and the space drag effect, is reminiscent of the Penrose effect. Indeed, the energy flux carried by the GM-field is 
nothing but the flux of angular momentum times lo, which completes the usual energy flux measured locally, as seen 
previously in subsection IjHI A|) . Infalling matter (n < 0) carries a negative specific angular momentum (Z* < 0), 
which, with the frame dragging, produces an outward energy flux. 

Following Michel , it is convenient to define the " magnetisation parameter" a as the ratio of the Poynting flux 
over matter energy flux: 

(35) 



awTdM 3 c 2 wTkc 2 

This number must be large at the beginning of the jet in order that the Poynting flux propells matter efficiently 
and then it diminishes in favour of the matter energy flux. The magnetization takes on a simple form at the Alfven 
surface: 

n»(fl* - aQ,-u)r A 
<?a = y~2 \ A > ( 36 

Far enough from the Black Hole, one gets a a — Qlr A /c 2 , and an Alfven surface far enough from the light cylinder 
would be expected, otherwise, the jet acceleration would almost end there. 

From the divergence free condition when radiation and dissipation can be neglected, a Bernoulli invariant can be 
derived (the fourth of the family): 

awTc 2 - -^-fi* + wluI ee B(a) . (37) 

Starting with Tq = (1 — ti^r^/ c?)^ 1 / 2 , B(a) = aQWoT (l + a a ) + wqujIq, with an initial magnetisation cr which is 
expected to be large in order that the magnetic field drive a relativistic flow. It is often convenient to choose these 
"initial" values at the light cylinder beyond the "watershed" surface. From (|37|l . it follows that, at remote distance 
from the Black Hole, the product T(l + a) is almost constant on each magnetic surface; and since Ta = k\I\ (where k 
is a constant), the Bernoulli invariant implies that r + fc|/| is almost constant on each magnetic surface. The expected 
growth of the bulk Lorentz factor requires a strong decrease of the electric current intensity, which cannot be described 
in force free MHD where the current intensity is constant on each magnetic surface. 

The decrease of the magnetization by a decrease of the current and an increase of the motions can be translated 
in angular momentum language, because of the ideal MHD assumption. Indeed the evolution of the magnetization 
parameter is described in term of the evolution of the angular momentum, by inserting the angular momentum 
equation IL'-'il) : 

- wl) 

° = a 2 Tc 2 ' (38) 
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the large initial value corresponds to a high value of the magnetic breaking effect with i < ( t , and the decrease 

at large distance where w ~ 1 corresponds to an increase of the matter angular momentum in widening magnetic 
surfaces that tends to /*. 

Far from the Black Hole, the Bernoulli function takes on a simple form, even after inserting the expressions for b, 
eq: and I, ea: pl : 

B(«) = wTc 2 - '»"» ~" m * f2 . (39) 
m z — 1 

the bulk Lorentz factor is then derived for prescribed magnetic surfaces : 

(m 2 - 1 + £, 2 )wT = (a + w )m 2 - w , (40) 

where £ = fi*f/c. We will discuss the asymptotic evolution of the Lorentz factor in the next subsection. 
The location of the Alfven surface is constrained by setting m 2 = 1, which leads to 

-^w a Ta = a . (41) 

For a not very hot plasma f& — t^{u o/F^) 1 / 2 ; at this stage however the value of Y a is unknown. Most of the Poynting 
flux would have been converted into matter energy flux at the Alfven surface and thus acceleration of the flow almost 
end there... Therefore the flow would be relativistic only if the Alfven velocity is relativistic, which amounts to the 
condition of a very dominating magnetic pressure there: B 2 /8ir 3> e + P. 

The non-relativistic limit (a ~ 1 4- <&/c 2 and w ~ 1 + h/c 2 ) of this invariant leads to 

1 o btt* 

-u 2 + <P + h - = B nr (a) , (42) 

Z 47TK 

where the mass contribution has been withdrawn. This is not necessary from a formal point of view, but convenient 
for the understanding of the difference between the relativistic and the non-relativistic regimes. Indeed this governs 
the motions along the poloidal magnetic field lines and particularly the increase of kinetic energy, which is easier 
than the increase of mass-energy in the relativistic regime for the Poynting flux. It is noteworthy to remark that a 
magnetisation parameter can also be defined in the non-relativistic MHD where energy mass is disregarded. That 
non-relativistic magnetisation can be made more easily large than in the relativistic case, which makes non-relativistic 
flow easier to be accelerated, moreover they do not suffer any radiation drag... But similarly to the relativistic regime, 
the acceleration of the flow is almost ended at the Alfven surface where tr^ r = 2J7 2 f 2 1 /V^|. Indeed Va/Q* is the typical 
scale of the MHD flow and the Alfven radius takes on values of this order. 

A fifth invariant is the plasma entropy (Pp~ 7 = K(a)). To close the system of equations that entirely describes the 
jet formation, an equation that determines the magnetic surfaces by stating the transverse balance is needed. This 
is the Grad-Shafranov equation (|3dj| These six equations determine a, &, I, u, p, P; the five invariants must 

be specified at some surface crossing all the relevant magnetic surfaces, e.g. at the Alfven surface or the surface of 
the rotor, and two boundary conditions must be set for the integration of Grad-Shafranov equation. Once a jet is 
launched, it has to cross three critical surfaces (see [73| and [8(| for a good analysis of the critical surfaces). In the 
case of the pulsar wind, the matter is so tiny compared to the intensity of the magnetic field, that only the light 
cylinder at = fi*/c requires a regularity condition. In AGNs, micro-Quasars and GRBs, matter pressure is not 
negligible compared to magnetic pressure. The Alfven surface replaces the light cylinder and there are, in addition, 
slow and fast magnetosonic surfaces. 

D. Asymptotic evolution of the flow 

From Bernoulli equation, the asymptotic evolution of the Lorentz factor of the flow can easily be derived; it is 
governed by the following equation when the flow has a non-relativistic temperature: 



T 3 -(l + a -a c +'pT 2 + ^- = 0. (43) 



In ea. l|43[) . the parameter a c corresponds to a kind of modified value of the magnetisation parameter that depends on 
the local poloidal field and is defined by 

q2 d r 2 „: 2 n 

(7 c = *7 . (44) 
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This parameter is close to the value of the magnetisation when the asymptotic magnetic configuration is described by 
a monopolar field: 

_ fllB r r 2 sin 2 9 

O~mono — . ; V 

ilTWKC 

in which B r r 2 = contant, thus leading to a constant magnetisation on the magnetic surfaces (cones in this approxi- 
mation) like in |4q . Clearly, a monopole approximation that would have been assumed too early in the description of 
the flow (this is an exact solution of the complete set of equations (^3|> 0) would prevent any significant acceleration 
since <r c ~ ao . Moreover an asymptotic monopole solution cannot exist if a c — ao as can easily be seen on eq. ( , 
that would have no root for f — > oo. Thus the difference Ac = oq — a c must be carefully treated in the flow description. 

This important equation ea.f l43f) has been obtained by inserting m 2 = u 2 /u 2 A — (u/a c )t; 2 and T 2 = 1 + u 2 + l 2 /c 2 f 2 . 
For large T, it turns out that the approximation u ~ T — 1/2T is sufficient. 

The fast magneto-sonic surface is another important critical surface located beyond the Alfven surface, but usually 
considered in the literature as not much farer. It can be found from ea.f I43[) which indicates that the derivative of T 
along the magnetic line has a denominator that vanishes when r = Tp with £ = £p such that 

3T 2 F = 2(1 + A<tf + <t c /Zf) ■ (46) 

The Lorentz factor at this surface must always takes the value Tp = a^ 3 (after having inserted ea. 1461) into ea. H43|) 
and the radius is such that 

(47) 



* L 3 1/3 1 A ' 

which shows that the FM-surface exists only if Act < \a\l 3 — 1. Therefore fp ~ f^Og 3 . Since Ta < Tp, it turns out 

1/3 

that ta ~ r Lcr also; these two estimates are general results for a cold flow. 
At this stage, a preliminary discussion can be addressed. 



If Act ~ ctq^ 3 even far beyond the FM-surface, then r m ~ a^ 3 and most of the Poynting flux has not been 
converted into jet power. The asymptotic solution can be of monopolar type. 

If AtToo ~ <t (which would necessarily held far beyond the Alfven surface only), a much better conversion would 
occur and Too ~ a®. Estimating ao at the light cylinder, the ratio between both magnetisations is 

ao _ \b^\ l ?1 = b ± B p r 2 | L _ n. - n L B p f 2 \ L 

a c TlB p f 2 TB p Bpf 2 B p f 2 

This ratio clearly tells how an efficient conversion of the Poynting flux into jet power could occur. The high 
intensity of the toroidal field compared to the poloidal one during the non-relativistic stage of acceleration is 
favourable, especially if the poloidal velocity u is still subrelativistic at the light cylinder (ul < 1)- But this 
is spoiled by the condition for smooth crossing the FM-surface (however see next paragraph where the case of 
a warm plasma is discussed) . The other favourable circumstance is based on the possibility that B 2 r 2 sin 2 9 
becomes sufficiently smaller than its value at the light cylinder. This is the attempt followed by N. Vlahakis 
(0). 

If the relativistic flow is collimated by a non-relativistic jet, such that > f^a^ 3 , and r 2 B 2 sin 2 9 < r 2 B 2 |j, 
sin 2 9 h-, then Acqo ~ o§. Therefore the flow undergoes a sizeable conversion of the Poynting flux into jet power. 

In the case of a warm flow. 



Crossing the FM-surface seems quite constraining for a flow having a non-relativistic temperature. From eq:(??), it 
is not difficult to modify the evolution equation of the bulk Lorentz factor for a flow that has a relativistic temperature, 
so that w = 47/3 > 1. It becomes 

wT 3 - (w + 0-0- a c w{l - l))r 2 + = . (49) 

One obtains the modified value of the Lorentz factor at the FM-surface Tp — ( ao+Wo ) 1 / 3 and the location of the 
FM-surface is such that 
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A relativistic temperature at the FM-surface greatly helps its smooth crossing, even for wo <C (Jo, provided that 
wp > 1 so that wpa c > gq. Then the FM-surface is located at a few light cylinder radii with 

fl=fl ^I^_ (51) 

WF&c - CO 

But the flow must have cooled at remote distance such that ~ 1 and Aa^ ~ ao to allow a large conversion 
such that Too ~ ctq. Thus the recipe to get a successful conversion of power is an initial magnetisation cr significantly 
larger than its monopole estimate, thanks to the non-relativistic stage of acceleration, together with a relativistic warm 
plasma up to the FM-surface. But the hotter the plasma the weaker the magnetisation; indeed for a relativistic plasma 
in mildly relativistic motions at the light cylinder, <jq ~ B 2 /8irP. In particular at equipartition the magnetisation 
parameter is of order unity and thus the formation of a relativistic jet becomes difficult. However a thermal relativistic 
wind is set up; this is the assumption of the "fire ball" model of GRBs (|z2) despite the very high intensity of the 
equipartition magnetic field that can reach 10 15 G... Nevertheless there also exist some theoretical att emp ts to explain 
GRB relativistic wind with the electromagnetic formalism presented in this section (as for instance [85j). 

It is quite conceivable that a cold flow with a high initial magnetisation can be entrained by a large Poynting flux 
without succeeding to establish a stationary flow that crosses the critical surfaces. This would give rise to Poynting 
flux dominated flow with a very large Lorentz factor for the electromagnetic field r e . m . 3> 1 with some amount of 
perturbing matter having T < T e . m . (Lyutikov & Blandford |4l|) • 



E. Blandford-Znajek mechanism 

The celebrated Blandford-Znajek mechanism (|Il|) consists in taping the rotation energy of a Kerr Black Hole. 
From thermodynamical arguments, the energy-mass of the black hole is shared in an irreducible part related to the 
entropy and a rotational part that can be extracted. At maximum rotation (a s — 1), Mj rr = and M rot — 

(1 — -^)M ~ 0.29 M. As seen previously, the energy extraction is done by the gravito-magnetic field ("space drag"), 
through a kind of Penrose effect. In the neighbourhood of the horizon, the effect of the magnetic field is supposed 
to dominate the inertial effect and a force free approximation can be assumed. The poloidal current is thus aligned 
with the poloidal magnetic field and therefore the current intensity is conserved on each magnetic surface: 1(a). Since 
the electric field is known, close to the horizon, the current is determined by the converging wave condition in the 
vacuum, as stated in subsection (3.2), namely B trans = —n x E(l + 0(a)), or — —Eg and Be = 0. From i(2*T)l . 
one straightforwardly gets 



ac 



f 2 B n (52) 



In agreement with the statements of subsection (3.2), both the transverse electric field and the toroidal magnetic 
field diverge like a -1 and contribute to the energy flux at 0(1). This form of b compared to (|20|l indicates that 
there is necessarly a constant flux of falling matter on the polar cap of the Black Hole such that apT = —n~B n 
(k — kT < 0) and v p = — 1 + 0(a). The total outgoing energy flux in the vicinity of the horizon is thus: 

P BZ = / a 2 S.nd 2 a = f QJQ H - ^*)r 2 — d 2 a . (53) 
Jn Jn 47rc 

This corresponds to the loss of Black Hole energy-mass per unit time in the form of rotation energy. This power is 
maximum when f2* = f2#/2, which can be interpreted as an impedance optimisation f |79j1. The electric circuit is 
supposed to be loaded far from the black hole by a dissipative region with a resistance equal to the resistance of the 
black hole, i.e. vacuum resistance. Some criticism of this optimisation has been addressed by Punsly and Coroniti 
r [&5| [66]), who argued that there is no causality relation between the flow beyond the matter source around the 
watershed surface and the flow falling toward the horizon and the ergosphere. However numerical simulations f |17| 
[l8|) have shown the capability of the process to generate relativistic jets, as long as no coupling with a radiation field 
is taken into account. 

As shown previously, the Poynting flux converges towards the Black Hole: as usual the Poynting flux converges 
towards the ohmic dissipation region, where the electric field and the current point towards the same direction. The 
dissipation of the electromagnetic field in the Black Hole is given by: 

P dls = [ (Q H - n,f . (54) 

Jn 47rc 
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This corresponds to the increase of Black Hole entropy as previously introduced (in subsection (3.2)), i.e. Pdis — 
Th<1Sh / dt. At the "impedance optimisation", Pbz = Pdis- At a few gravitational radii, the space dragging decreases 
enough so that u> < O*, the electric field changes its direction, and the Poynting flux dominates, carrying energy- 
momentum far away and eventually gives it up to matter. This is potentially able to generate a relativistic jet of 
power: 



It is interesting to compare it to Eddington luminosity and also to bear in mind the Eddington scaling of the magnetic 
field: 

Pbz ~ m-4 2/ B \2/ M . 4 2 . P , 2 f r a \ 

~L~Ed as W ( To^ } " 10 aA p7? ■ (56) 

This ratio is independent of the Black Hole mass and it is not expected that the magnetic field could be larger than ten 
times -B*, the Eddington estimate for a s = 0. It clearly shows that the Blandford-Znajek mechanism has a relatively 
significant power only for a central engine that works well below the Eddington luminosity, as would be the case of a 
weakly radiating torus. 



V. RELATIVISTIC AND NON-RELATIVISTIC JETS 



How a Poynting flux is converted into energy-momentum of matter? This is a difficult problem that has not 
been solved yet, although there are some possibilities previously discussed. This flux propagates through a plasma 
which, first, is falling (k < 0) on the Black Hole, then is generated in the form of a pair plasma around a stagnation 
surface, and finally expelled (n > 0) (see figure (1)). Looking at the Bernoulli invariant (|37|) . it can be seen that the 
asymptotic bulk Lorentz factor is directly determined by the value of the Bernoulli function, especially if the current 
decreases along the jet, which is likely. As previously seen, to get a high Lorentz factor, it is necessary to have a large 
magnetisation, which means a strong contribution of the magnetic field to the angular momentum transfer, at the 
beginning of the jet since, 

13(a) = auTc 2 (l + a) + wojl (57) 

Moreover the strength of the magnetisation is nothing but the strength of the current (or toroidal field) at the beginning 
of the jet, which constrains the angular momentum invariant. It is thus expected that there exists a starting region of 
the jet fairly free from gravitation, such that the total angular momentum is dominated by the field, i.e. |&|/47tk ~ /*, 
and (j — <t* = Z*f2*/c 2 . Thus the conversion of the Poynting flux would be successful and c 2 Poo ~ B(a) with B ~ 
This is actually its maximum possible value; however no complete analytical nor numerical solutions have been derived 
so far displaying such achievement. This is probably due to the defect of the "monopole" type solution that has been 
used as a starting solution, because it is an exact solution of the complete set of equations (as explained in subsection 
(4.4) and also N. Vlahakis's chapter). 

Similarly, in the non-relativistic case, the asymptotic bulk velocity is approximately given by the value of the 
Bernoulli function: ~ \J2B nr . Similarly to the relativistic case, the Bernoulli function is determined close to the 
source where the magnetic field is expected to dominate: B ~ — £!*j^t|o ~ f2*Z*. Inserting the angular momentum 
equation i|23|) in the non-relativistic Bernouilli invariant 142|) . one gets : 

l-u 2 + $ + h + n*(l*-l) = B nr (a) , (58) 

Since / — *■ Z», the same result is obtained for both relativistic and non-relativistic cases: B nr ~ Z»fi*. In both 
case, the matter angular momentum starts with a low value compared to the toroidal field contribution and then 
increases whereas the toroidal field contribution decreases. Therefore the relativistic flow is achieved if Z*J7* > c 2 and 
Too ~ Z*r2»/c 2 ; if not, a non-relativistic flow is obtained with Uoo ~ \/2Z»f2*. This seems simple, but the way these 
invariants Z» and are determined is not simple. 



A. Compton Drag 



The main danger for the relativistic flow is the Compton drag. It is easy to estimate the characteristic length of 
relaxation of a cold relativistic flow that suffer energy-momentum loss by Compton radiation. Indeed, assume a mass 
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flux Mj of the relativistic flow, composed of electrons, positrons and protons satisfying electric neutrality, with an 
electron density n e = n (comoving) , a positron density n + — (1 — £)n and a proton density n p = £n. The mass flux 
across dA is 

dMj = ((2 - £)m e + ^m p )nT(3cdA 
The energy decrease due to Compton drag can be described by the following equation: 



div{S m + S e - m - g -) = --(2 - OnvTcU^T 2 < 7 2 > , (59) 



where 7 is the Lorentz factor of an electron measured in the co-moving frame and U-y is energy density of the 
radiation field. This expression of the Compton drag is obtained with the approximation of an isotropic radiation 
field; anisotropy effects will be analysed in a specific subsection (7.1) devoted to the Compton rocket effect, which 
necessarily requires a kinetic description. A relaxation length £ re i(r) can be properly defined such that 

div[{l + a)S m ] = -(1 + a)S m /£ re i . (60) 

One finds 

Uei{r) = 1 y a -rz -z C + 5 — r , (61) 

where the compactness £ c , as introduced in first section, is the main parameter that must be compared with the 
magnetisation a. The relativistic jet is killed by the Compton drag if £ re i < r. Two important cases must be 
considered: the electron-positron jet and the electron-proton jet. 
In the case of a pair jet, 

1 + cr < 7 > 

4ei = —„ * — r . 62) 

4 < 7 2 > 

as can be seen, a hot distribution of pairs makes the Compton drag more efficient than a cold distribution. The 
relativistic jet can be launched only if <tq ^> 4(^o) < 7 2 > / < 7 >■ This constraint is severe and usually not achieved. 
Indeed, as previously seen for a luminosity smaller but still of the order of Eddington luminosity, the compactness 
4 = 10 2 - 10 3 . It turns out that pair creation can be catastrophic with high compactness and unavoidable in situ 
reacceleration of them. Therefore they get a pressure that easily reaches the magnetic pressure, which makes a ~ 1. 

The relativistic jet immediately relaxes by emitting hard X-rays. 

1/3 

Note that there is an amazing regime when <jq 1 and a < 4 < Co, that has been pointed out by Begelman 
et al. ([13) and Beskin et al. (0|), because of a non-trivial effect of the magnetic field line curvature. In this 
regime, although the Compton drag reduces the power of the jet, the bulk Lorentz factor is increased compared to 
the Compton drag free case. 

It is often thought that a large proton load of the flow could help, but even if we assume that there are as many 
protons than electrons at the origin of the jet, the relaxation length is such that: 

^ - i±^!^| r . (63) 

4 m e < Y > 

The mass ratio makes it more favourable, but still a hot distribution of electrons makes launching difficult. Moreover 
protons need to be accelerated before competing with the Compton drag. In the next subsection, the proton entrain- 
ment will be estimated. Anyway the important conclusion that can be drawn about the possibility of jet production 
by the Blandford-Znajek mechanism is that it requires a low compactness, which is possible only if the accretion flow 
is not radiatively efficient (L < 10~ 3 LEdd)- 



B. Baryonic pollution of a pair jet 

Diffusion of protons from the accretion disk or torus or from the non-relativistic outflow towards the pair jet takes a 
long time, even if we consider the fastest diffusion coefficient across the magnetic surfaces, namely the Bohm coefficient 
(Dj_ = ^vtl). Remark that the thermal estimate of Bohm diffusion is the same for both electrons and protons if 
they have the same temperature (cT/eB). Assume a constant flux of baryonic matter crossing the magnetic layer of 
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radius a and length Ar that interacts with the leptonic flow that fastly entrains it. Then in permanent regime, the 
ratio of the number density of polluting protons over the external protons is given by 

n mt 2D ± Ar t translt 

(64) 



n ext Tv p a 2 tdijfus 

where v p — [3 p c is the poloidal velocity of the jet. This can be rewritten in term of the travel length Ar m i X required 
to get full diffusion, n mt — n ext : 

. 3 Tc a 

Ar m « = -— a , (65) 

2 p p V A r 

where tq = Va/oj cp ; that kind of typical Larmor radius is the minimum scale of collisionless MHD description. It 
turns out that a/ro is a huge number (which, by the way, validates collisionless MHD in these topics), indeed 

T JL * 0.6 x iq" (_5L_^/a . (66) 

ro lO a M0 10 1/£ cm 6 

The density ratio is thus 

n ext Ar mlx 

For instance at 10 3 r 9 , this ratio is less than 10 -10 (the ambient medium has a density much smaller than 10 12 cm~ 3 , 
10 9 cm -3 say), the entrained proton density is ~ 10 _1 cm -3 , whereas the pair density is ~ 4 x 10 3 cm -3 (when the 
pair pressure is in equipartition with a magnetic field of 1G and an average energy jm e c 2 with 7 = 10 3 ). 

This issue of baryonic pollution of a relativistic pair flow is particularly important in the physics of GRBs also. 
In the fire ball model, the asymptotic bulk Lorentz factor is determined by the amount of baryons that have been 
entrained by the e + — e~- wind; and the model suffers from requiring a fine tuned value of the baryon parameter to 
account for the observations, without any physical argument. 

Note that the pair flow is first polluted by the most energetic protons that have a larger diffusion coefficient through 
their larger Larmor radii. 

Another way to accelerate a relativistic pair jet will be presented in a next section devoted to the "two-flow" model, 
where the Compton rocket effect in a hot pair plasma is taken into account. 



C. Confinement 



Because of the toroidal field, it may be expected that the flow tends to self-collimate, if the magnetic tension 
compensates the outward pressure force and the centrifugal force. A necessary condition for self-confinement is that 
the total value < 0, which reads: 

(e + P)T 2 n 2 f 4 < — (68) 

47T 

The toroidal field, and thus the current, decreases with mass density along the jet. Because I ~ I* for m 2 3> 1, the 
collimation requirement becomes asymptotically 

6 2 /4tt > pll (69) 

Although self-collimation can be realized in non-relativistic jets, numerical works have shown a lack of collimation of 
relativistic jets (01 01)- This suggests that the tension is unable to balance the internal pressure in relativistic jets. 
The difference between relativistic jets and non-relativistic ones comes from the importance of the electric field in the 
former case. Indeed B 2 , — B 2 (Q„ — Q) 2 f 2 /v 2 and E 2 = B 2 £l 2 f 2 /c 2 . This clearly shows that, for v p — c, the toroidal 
magnetic field is a little smaller than the electric field, and, thus, its tension is unable to balance the pressure due to 
the electric field. 
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On the contrary, with non-relativistic motions, cylindrical self-collimation is possible, as rigorously proved by Heyvaerts 
and Norman [33. Even recollimation is possible after the flow has passed the fast magnetosonic surface (\58\. Hi|), 
thus assuming a tulip shape. 

These statements can be set more accurately. Suppose that a cylindrical collimation is realised with a necessarly 

finite amount of current intensity 1(f) — > I to t at large distance from axis. Let P' = P + — £ g^— -; then the confinement 
is such that the toroidal field contribution balances the gradient of P' and the centrifugal effect: 

% -^Wf-~±4f<rB^. in) 

First, suppose that E 2 <C P 2 , as it is the case in non-relativistic jets; then one derives from the confinement equation 
the required current: 

L^l= [ (2rP' + wpT 2 Q 2 r 3 )dr-f 2 P'(f) (72) 
2tt Jo 

B 2 , P' and p are supposed to radially decrease so that the integral converges, and thus the existence of a jet radius 
means that at larger radii the current saturates at a finite value. 



-^L = / (2rP' + w P r 2 n 2 r 3 )dr . 
2n Jo 



(73) 



A jet radius fj can therefore be defined such that = 47rf|P'(0). 

Now, consider a relativistic jet developing an electric field such that E 2 > B 2 ^. Take the estimate E 2 = B 2 ^ + 2B 2 
(in fact, the poloidal contribution is not important, as will be seen below). Then the confinement equation is changed 
into: 

f)P" 1 Ft T 2 

~ w P T 2 n 2 f = -^(p 2 - U • (74) 
or irr z or r 

The modified pressure term P" now contains 3B 2 /8tt instead of B 2 /8tt. The solution is such that 

L^l = f [ (P" + wpT 2 n 2 r 2 )dr - f 2 P"(f) . (75) 
Jo 

Since collimation is assumed, the integral must converge and f 2 P"(f) — > 0. But then / should be proportional to y/f, 
which is in contradiction with a finite confining current, since the existence of jet radius fj would be characterized by 
a convergence of the current for f > fj. Therefore it is impossible for a relativistic jet to undergo self-collimation. 
As anticipated above, the modification of the pressure term (P" instead of P') does not make a qualitative difference 
and, even with a lower estimate of the electric field, E 2 = B^, collimation is impossible. 

To achieve self-collimation, the jet must be composed of a relativistic flow (hereafter R-jet) close to the axis 
embedded by a non-relativistic one (hereafter NR-jet) that insures collimation. The relativistic is confined in a radius 
r* and the non-relativistic one in a radius fj and they are related to the currents through the following relations: 



and 



I 2 (r* 



2tt 2tt 



(P" + w P T 2 tt 2 r 2 )dr - riP" (r*) , (76) 



I (2rP' + w P n 2 r 3 )dr + r 2 P' 

J r* 



(r.) . (77) 



By inserting profiles for the pressure and the mass density in these integrals, one gets a relation between the radius 
and the current. For smooth profiles in the relativistic jet, one roughly gets the relation 

^(e*+P^ 2 (r*)-P>*) • (78) 

And for profiles like ( I \ ) 2 for P' and p in the non-relativistic jet, one gets, for rj ^> r* 

-f-r 

I 2 ot ~ 7 2 (n) + 2n(P'(n) + \p(r»)ttt{r*)rl)r 2 . (79) 
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Since the angular momentum is related to the initial current through Z* ~ —Iq/2ttk, and reminding that Attk 2 = 
Pa/wa, one has 

^ - ^< l* • « 

where <S» labels the interface between the R-jet and the NR-jet. This stresses one more time the necessary decrease 
of the electric current in jet generation. Another important remark is that the power of the NR-jet needs to be 
significant in order to confine the R-jet, since the power is proportional to the current, at least in the launching stage. 

VI. JETS FROM AN ACCRETION DISK 

The accretion power is closely related to the process of angular momentum extraction from the disk. Powerful jets 
can be generated from the accretion disk provided that they play a role in the angular momentum transfer. There 
are four possibilities for transfering angular momentum of a disk. (In this section, the cylindrical coordinates (r, <p, z) 
are used) 

• The usual viscous transfer: 

an 

which is very unefncient to account for the observed disk luminosity. 

• The transfer by the MHD turbulence stress, as generated by MRI for instance: 

• Instead of being generated by an MHD turbulence, this stress can be generated by a large amplitude wave with 
m = 1 excited by the spiral wave instability when the magnetic field is above the equipartition value ([j^ 

• Instead of transfering radially, a static magnetic field with opened lines in a bipolar configuration exerts a stress 
that carries the angular momentum vertically: 

This last possibility deserves a specific investigation, for this is precisely what jets are doing. 

A. Angular momentum transfer by the jet 

It is straightforward to compare the torque due to the static magnetic field, that has both a significant vertical 
component and toroidal component, to the one exerted by MHD turbulence in its most efficient regime where r^ r ~ P m 
(see the section devoted to the magneto-rotational instability). The torque per unit volume is r'S/jT^. Thus the ratio 
of the turbulent torque over the static field torque is 

turb torque hB~ 

(81) 



field torque rB^ 

It is easy to get \B<j,\ > —B z . Thus the vertical extraction of angular momentum by the jets is easy compared to the 
radial transfer by turbulence (|H3| Actually, the main problem will be to avoid a too strong generation of B^ 

by differential rotation in the accretion disk by an appropriate resistivity (this is likely an effective resistivity that 
accounts for reconnections in the disk) (|2{|). 

The transfer of angular momentum relates the accretion rate to the torque (bearing in mind that b(r, —z) = —b{r, z)): 

M a (l -li) = - [ b + B z rdr - 2irr 2 [ r^dz (82) 
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The Shakura-Sunyaev prescription consists in assuming r^ r = —a v pil 2 h 2 or a viscosity of the form v v — a v flh 2 , 
and it is expected that the turbulence allows to reach a value of a v as high as 0.1. Similarly a turbulent magnetic 
diffusivity v. m = r\j^K can be produced such that v. m — a m Qh 2 , with a m ~ a v . The angular momentum extraction 
by the static field can be characterized by a parameter fi such that B^B z /A-k = —ppQ 2 h 2 ; \i is supposed smaller than 
unity, as well, for a reason that will be explained later on. The accretion rate that stems from both kinds of angular 
momentum transfer is thus 

M a ~ 4ir(a v h + pr)pVth 2 . (83) 

A magnetic Reynolds number lZ m can be defined as lZ m = ru^/v m , where uq is the accretion velocity and its value 
is lZ m ~ a ^ e+ J 1 where e is an appropriate average of h/r. 

Assuming that the magnetic stress extracts a fraction \j 01 the angular momentum and the turbulent stress a 
fraction Xd {Xj + Xd = 1), then 

• 81 

rb+B z = -X 3 M a - (84) 



whereas 



^-(2^r 2 j T<t>r dz) = ~XdM a ^- (85) 



The power budget is therefore, as follows; for each jet: 

Xj 1 GMM a 



Pi = (86) 

2 2 r t K J 

and for radiation: 

1 GMM a , N 

L = Xd- 2 —^ (87) 

In fact the accretion power is shared in several contributions, that are simply expressed in terms of the Reynolds 
number lZ m and the magnetic stress parameter \i. When the magnetic Prantdl number is supposed to be 1, for 
simplicity, we have : 



a Joule heating contribution 



a viscous heating contribution 



therefore a luminosity 



• and a power for each jet 



where the Reynolds number 



Pjoule ~ ^-Pacc , (88) 



T^Pacc , (89) 



^acc j (90) 



K 



Pjet-l(l-^-)Pacc, (91) 



n m ~ i + ^ . (92) 

e 



Thus disregarding the small contribution of the Joule heating, one has \d ~ and \j ~ 1 — 
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There is a relation between the accretion rate and the ejection rate derived from the conservation of the angular 
momentum flux in the jet, indeed the angular momentum flux in the jet, l^dMj , comes from the fraction of the fraction 
of angular momentum extracted from the disk by the magnetic stress : 

dMj = |M a ^ ; 

which shows that the widening of the magnetic surfaces allows to transfer the angular momentum with a rather weak 
ejection rate, provided that the Alfven radius (the "lever arm") be large enough, since, roughly, Mj/M a ~ Xj r \l r \- 

The first solutions derived by Fcrreira & Pcllcticr ( 25] ) were obtained even with the extreme assumption of no 
other angular momentum transfer than through the magnetic stress, which corresponds to > e. The conditions 
for this accretion-ejection process to work are quite stringent but sensible. They have been made more flexible by 
introducing also an effective viscosity and an additional heating at the surface of the disk j2jj) and checked by 
direct numerical simulations (|2lj|). It turns out that without this additional heating, \d ^ lj most of the accretion 
power goes in the jet, whereas the additional heating at the surface makes the mass outflux easier and allows solution 
with sizeable conversion coefficients 0.1 < Xd ~ Xj < 1- 

Such a jet that takes its power from extracting a significant fraction of the disk angular momentum has specific 
properties that have been investigated in( [53]) and allow observational tests (|63|)- 



B. Power and asymptotic motion 

There is a natural scaling with the mass of the central Black Hole for the power and the mass flux of the jets 
launched by accretion disk. The mass ejection rate is a fraction £ of the accretion rate Mj = £M a ; this ratio depends 
on the details of the transport physics in the accretion disk but not of the Black Hole mass. The averaged asymptotic 
bulk motion is such that 

r-i = ^^ = ^^_ (93) 

M jC 2 ^n so 

The bulk motion is independent (or weakly dependent) of the Black Hole mass. This is probably the reason why 
similar motions are observed in quasars and micro-quasars. This argument favors the model of jet generation by 
the accretion disk. However, because r — 1 ~ 2 x lO -2 ^ -1 , £ ~ 10 -3 is required to get a relativistic motion with 
r = 10; which is difficult to obtain from an accretion disk (|25j); higher values are commonly obtained in numerical 
simulations. However relativistic motions are not ruled out in the innermost region and could be obtained from a 
weakly radiative torus. These estimates must be linked with the results obtained with Bernoulli equation (see section 
0)- This relates the angular momentum invariant with the ration £ since — 1 ~ ~ 2 x 10 -2 £ -1 . 



C. Conditions for ejection 

The accretion flow bends the field lines and the generation of a strong toroidal field by differential rotation must 
be prevented by the appropriate resistivity. Assume a resistivity profile v m = u m o(l ~ i?) with a kind of Sakura- 
Sunyaev prescription u m o — a m flh 2 , then the vertical variations of the field components can be expanded in Legendre 
polynomials. In this idealised description, the matching between resistive MHD inside the disk and ideal MHD in the 
jet has to be set at the disk edge z = h(r). The value of a m directly controls the relation between toroidal and radial 
components at the surface of the disk, B~f and B+ : 



The parameter A is a pure number that depends on the geometry of the field lines inside the disk, which is the 
solution of a kind of nonlinear eigen- value problem ( 57]). For instance, the less "wavy" solution, where B^ varies 
almost linearly with z across the disk, leads to the lowest value of A (ground state), which is close to unity. However, 
because a m is expected to be significantly smaller than unity and not significantly stronger than the other field 
component, in order to prevent a strong squeezing of the disk, a higher value of A ~ 10 is better, which is obtained by 
the second eigen-state, where has a cubic like variation across the disk. The poloidal field is quartic rather than 
quadratic in the vicinity of the equatorial plane. 
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It is expected that the field lines bend above the disk in order to fulfill the Blandford & Payne (H2]) condition for 
launching in ideal MHD regime, namely a bending of more than 30°. Moreover it is expected that the field does not 
squeeze the disk significantly: P m < P. In fact, it is possible to get B~£ ~ _B+ ~ £>+ together with a m ~ 0.1 

In the neighbourhood of the disk surface, matter is expelled provided that: 

d$ eff dP 1 db 2 n 
os os 8irr z os 

where s is a coordinate along the poloidal field line that vanishes at the foot of the line. Expansion of the effective 
potential for small s leads to — ^o(3 tan 2 ip — l)zcosip, where ip is the angle of the tangent to the field line 

with respect to the vertical direction. Therefore, if the field structure is almost force free (b(a)) and matter cold, the 
condition for an outflow is that the field line bends by an angle larger than tt/6. The effective force vanishes at the 
point where ip = 7r/6. Consider now, finite temperature and non force free effect. If at this point, MHD is already 
ideal, the current is decreasing and acts like the pressure gradient to expell matter. But if MHD is still resistive (inside 
the disk), the toroidal field squeezes the disk and tends to prevent matter escape. Thus the pressure gradient must 
overcome the squeezing effect; the condition is roug hly P > (B+ 2 + B+ 2 )/8tt or simply P > P m . This is the reason 
why there is a compromise for the intensity of the magnetic field: it must be high enough to extract the angular 
momentum, but the generation of the toroidal field must remain moderate to avoid the disk squeezing that prevent 
matter escape. This compromise is a magnetic field intensity of the order of equipartition value. Clearly the previous 
condition indicates that a hot corona helps matter to escape. When an additional coronal heating is at work ( |2C|p. 
the mass ejection rate is related to the coronal temperature. 

Although the first semi-analytical solutions (j2fi|) to the accretion-ejection problem were derived with a self-similar 
assumption, that cannot be valid at both edges of the accretion disk, they revealed the main constraints for an 
accretion disk to launch jets. The main features of these analysis were confirmed by direct numerical solutions (|2jj). 
Now to get fair simulations of the asymptotic behaviour of jets, 3D-simulations are needed to understand how the 
development of 3D-instabilities does not destroy the jet (H3, 0). 



D. Resistivity problem 



In a SAD, an effective viscosity is invoked to get the expected accretion rate to account for the luminosity of the 
disk. In AGNs for instance, the viscosity coefficient a v of the Shakura-Sunyaev prescription (e.g. T<j, r = a v pQ 2 h 2 ) 
is expected to be on the order of 0.1 and supposed to stem from a turbulent state of the accretion flow. If the 
angular momentum is vertically extracted from the disk by the magnetic field, as explained previously, then there is 
no viscosity problem. But there is a resistivity problem, because a quasi stationary accretion flow cannot be set up if 
matter does not easily cross the magnetic surfaces. A magnetic diffusivity is required as expressed by the prescription 
Vm = a m ^h 2 with a m ~ 0.1 as well. In fact, a fully developed state of MHD turbulence leads to turbulent magnetic 
diffusivity comparable to the turbulent viscosity. But as shown before, even if such a turbulent transport is achieved, 
the vertical transport of angular momentum by the jet can easily dominate the accretion dynamics. 

Knowing the cause of the excitation of turbulence in an accretion disk is an important issue that receives crucial 
contribution during the last decade. 



E. Magneto-Rotational Instability 



The magneto-rotational instability has been considered as a good driver of the turbulence in the disk, with the 
capability of transferring the angular momentum radially [HI ; 0)- A brief presentation of the Magneto- 

Rotational Instability (MM) will be done, using the formalism previously introduced in section (3) and subsection 
(4.1), and will be limited to the following assumptions for the sake of simplicity: the analysis is done in a thin 
equatorial disk far enough from the Kerr Black Hole where u> ~ 2a s cr 2 /r 3 and Hr^ ~ ui (see subsection IIII A|) . with 
a vertical mean field Bq(t), with incompressible perturbations. The main purpose is to look at whether both MRI- 
instability and jet launching are compatible in an accretion disk. The radial equilibrium for ZAMOs is such that the 
centrifugal force modified by the gravito-magnetic effect balances the gravitational attraction g r = — c 2 ^ log a with 
a 2 ~ 1 — 2r g /r; which reads 



fl 2 + ljQ = SI 2 , = -c 2 -— log a 
r Or 



(95) 
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Consider a small radial displacement £ such as au r = Since the magnetic surfaces are frozen in the flow as 
expressed by eq. , the linearised version reads : 

a = 5a + £— = 
Or 

The evolution of the radial displacement (Alfven wave modified by centrifugal force) is governed by: 

|^-V2|^-r(2fi + W )«l = 
The angular velocity undergoes Lagrangian variations 57 such that: 

n = sn + ^ . 

Or 

Note that the Lagrangian variation of the specific angular momentum I = 51 + differs from Or 2 by the Coriolis 
contribution 2Qr£. The equation governing the radial displacement can be rewritten as: 

8 2 8 2 Bfl 

This form of the equation already shows that the angular velocity shear is the motor of the instability. The coupling 
with the Lagrangian variation of the angular velocity will introduce more inertia due to Coriolis force. A second 
equation coupling Q with £ is obtained by the transport of angular momentum and the induction equation, by 
eliminating the toroidal field b between them: 

m = i^al (96) 

db gOr 2 2 dujd£, , 

The second equation describing the MRI instability is thus obtained: 

Looking for eigen-modes and defining the appropriate vertical wave numbers k n , which has a minimum value *~ tt/Ji, 
the following condition for the development of the MRI instability is derived: 

k 2 n V% + (2n + w)r-H-(n + u) < . 
or 

The condition can be rewritten in the form : 

klvl < 3^(1 + T JL + 2a s (^) 3 / 2 + ...) . (99) 
r r 

This result is close to Charles Gammie's result (32]), it differs by the sign of the Black Hole rotation contribution 
in a, which increases the instability in this result. The result is obtained by expanding in r g /r, however it seems 
that the instability is stronger at shorter distance. The interest of the exercise is two-fold. First, disregarding the 
gravito-magnetic contribution, it shows that there is a threshold condition for the growth of the instability in an 
accretion disk, because of the minimum value of the wave number. There is an instability only for B < B c such that 
B 2 /8ir ~ P. This condition is more or less in contradiction with the condition for jet launching. In other words, 
the interval of magnetic field intensities in which both the MRI develops turbulence and jets are launched is narrow. 
Second, the gravito-magnetic contribution in the vicinity of a Kerr Black Hole enlarges this interval and strengthens 
the instability together with the Schwarzschild contribution in r g /r. Anyway another instability would be welcome... 
A detailed analysis of the MRI instability in the vicinity of a Kerr Black Hole will be found in a forthcoming paper 

(H2). 

The MRI has a maximum growth rate of order Q, for k max — Q,/Va, if larger than n/2h, and excites 3D-turbulence 
at scale h and leads to a stress r^ r ~ P m . The resistivity is such that the effective magnetic diffusivity remains below 
but on the order of the value which would remove the frozen in condition for the mean field; 

Vm ~ n/k 2 rnax ~ v%/n ~ (p m /p)nh 2 , 

which is the expected height, but is valid for P m < P, and corresponds to a Shakura-Sunyaev prescription with 
a v ~ ct m - P m /P < 1. 
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VII. TWO-FLOW MODEL FOR AGNS AND MICRO-QUASARS 

The two-flow model consists in combining a relativistic flow along the axis and a non-relativistic or mildly relativistic 
flow coming from the accretion disk. When this model has been proposed in the eighties (mostly with a description 
of the interaction of the cold relativistic beam with an ambient medium) ([5(|), the arguments were rather poor: 
the phenomenology of hot spots do not require a relativistic flow, whereas relativistic sporadic ejections are observed 
only at short distance from the source; the large scale jet does not seem to loose significant energy during its travel, 
and non-relativistic shock acceleration satisfactorily explains the hot spot synchrotron spectra. So the large scale 
NR-jets and the short scale Rjets could be generated differently. First astrophysical applications of the model was 
published in 1989 (|75| |59| [60| ) - But more arguments have been proposed later ([3{|, ,36]) in relation with X-ray 
and gamma-ray emissions. The interest of the two-flow model relies on the following main arguments: 

i) The relativistic jet alone would radiatively cools rapidly and be relaxed by Compton drag (|S3). As we previously 
saw, Blandford & Znajek mechanism works with the bunch of magnetic surfaces that thread the Black Hole ergosphere 
thanks to falling matter carrying negative total angular momentum (mostly due to the toroidal magnetic field). That 
matter is necessarily created above the polar cap, and cannot be anything else than pairs. Above the stagnation 
surface the outflow is made of pairs. They are unavoidably experiencing the Compton drag. People claim that there 
would not be any Compton drag problem with protons. But protons can come in the relativistic jet only by diffusion 
from its environment (accretion disk or MHD NR-jet) and then be entrained. This amount of protons is very tiny up 
to a large distance beyond which ambient protons fill up the relativistic flow. 

ii) The relativistic jet does not self-collimate. The non-relativistic flow, which self-collimates, can confine it and also 
manages an ecological niche for it, because of the centrifugal barrier inherently associated with the MHD acceleration 
of the NR-jet. A simulation of a collimating shock between the R-jet and the NR-jet has been done ( 80]). 

iii) The extraction of rotation energy from a spinning Black Hole does not generate a power strong enough to 
account for the large power of FR2 jets. Jets powered by accretion are more powerful, but are mildly relativistic. 

iv) The R-flow is likely in the form of flaring clouds, rich in e + — e _ -pairs, channelled by the NR-jet, and involving 
a smaller power. It probably does not travel beyond kpc distance from its source. 

FIG. 2: The two-flow model (from L. Sauge thesis manuscript). The non-relativistic MHD jet, launched by the accretion 
disk, manages a centrifugal barrier along the axis where pair clouds can be generated by photon-photon interactions. In the 
anisotropic radiation field emanating from the accretion disk, these pair clouds suffer a Compton drag or propulsion depending 
on their velocity and energy distribution, and they are channeled by the magnetic surfaces. They are maintained with a hot 
relativistic distribution by in situ acceleration through the MHD jet turbulence. 



A new version of the two-flow model has been developed since 1991, where the R-jet is a hot pair plasma, energised 
by turbulence from the NR-jet (H^, [Ii|) and experience a Compton Rocket effect. Satisfactory fits of the spectra of 
3C273 and 3C279 has been obtained, with a natural break at a few MeV (0])! the stratification plays a major role in 
the analysis, which deeply differs from homogeneous or slab models. Later, a study of the micro-instabilities of such a 
R-jet have been investigated (0); these instabilities differ from the one investigated in the previous version, by the 
fact that they are triggered by the streaming of the cold ambient medium in the hot R-jet, instead of considering the 
R-jet as a cold beam exciting waves by pervading the ambient medium. This is an interesting circumstance of plasma 
turbulence excitation, because the streaming instability is driven by the radiation force and quasi-linear theory is 
not sufficient to describe the stationary state; non-linear mode couplings are necessarily taken into account, and, in 
particular, allow to estimate the heating of the pair plasma. 

A. Compton Rocket 

Following O'Dell (|55|). the anisotropic radiation field of soft photons, emanating from the accretion disk, or from 
the torus, generates a Compton force on electrons and positrons that can be pushing or dragging depending on the 
velocity of the flow. The component Fq = F.v of the force depends on the three Eddington momenta of the radiation 
field, on the velocity j3 of the radiating electron and its pitch angle cosine /i; it is such that: 

J = F = -2ttct t7 2 /3[/3(3J - K) - 2(1 + f3 2 )H^i + 0(3K - J)^} , (100) 
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where J, H and K are the Eddington coefficients, momenta of the intensity function : 

J= J I(v,n)dvdii s , 

H = ^ J J I(v,n)fi s dvdn s , 

K=^J J I{v,n)ii 2 s dvd^ s . 

These momenta have been calculated for a SAD and there exists an equilibrium velocity j3 eq that makes the bracket 

in F to vanish, such that 4| = (H^|, 0> H3)- For P < Peqt the radiation force accelerates, because the head 
on collisions of the electron (positrons) with photons favour a forward momentum in statistical balance, whereas 
it decelerates for (J > j3 eq (Compton drag), because electrons (positrons), advancing at relativistic motions, suffer 
more recoil effect in the statistical balance. For a SAD, T ~ {f-) 1 ^ 4 when r, <C z <C r e , and r e respectively 

being the internal and external radii of the SAD, and V ~ (— ) when z ^> r e . A detailed analysis of such a 
propulsion has been done by Rcnaud and Henri ( 69] ) for AGNs and micro-quasars. The results for micro-quasars 
cannot simply be obtained by a black hole mass re-scaling because the Klein-Nishina regime of the Compton effect 
is more important than in the case of AGNs (see fig(3)). Clearly, the Compton rocket is less efficient when the 
Klein-Nishina regime becomes significant. Moderate bulk Lorentz factors were obtained, but they are sufficient to 
account for the observations. In fact, a recent study by Henri and Sauge shows that the blazar bulk Lorentz factor 



FIG. 3: The jet bulk Lorentz factors (from 69]). The terminal bulk Lorentz is obtained as a function of j ma x, which is more 
or less the mean energy of the pairs in the co-moving frame, for different values of the distribution index s. The two top panels 
corresponds to a 10 9 M©-black hole with r e — 3 x 10 3 r 9 , L = LEdd (left) and L = Q.lLEdd (right). The two bottom panels 
correspond to a 5M0-black hole with r e = 3 x 10 3 r 9 , L = LEdd (left) and L = O.lLEdd (right). 



should not be very large, larger than 2 but significantly smaller than 10 ( 72]), contrarily to the estimate derived from 
homogeneous models. Indeed homogeneous models provides Doppler factors larger than 30 when the correction of the 
Infrared background absorption is done (even by lowering the doubtful point at 60/im. These high estimates are in 
disagreement with the reasonable assumption that BL-Lac objects are normal FR1 jets seen at narrow angle, which 
leads to statistical estimates of bulk Lorentz factors between 3 and 5. Stratified models where the X-ray synchrotron 
emission is largely due to other electrons than those emitting the inverse Compton emission can fit the compound 
spectra with these low values of the Lorentz factor. However they require local quasi mono-energetic distribution 
functions rather than power law distributions and an opacity implying a significant pair creation (see |6l|h These 
observational constraints revealing the importance of pair creation has a tremendous impact on the generation of the 
jets, as emphasized in several sections. The fact that the local distribution functions are likely quasi mono-energetic 
rather than of power law type implies that the acceleration process at work is more likely a kind of second order Fermi 
process or due to an electric field in reconnection sites. 



B. Variability 



The physics involved in the two-flow model makes the relativistic jet intrinsically variable. If the turbulence heating 
of the pair plasma is maintained in a region where r 77 > 1, pairs create 7-rays by inverse Compton effect, that are 
converted into new pairs, which are accelerated and then create also 7-rays, and so forth. This is the pair creation 
catastrophe, which leads to a highly unstable situation. Its simulation shows various regimes of flaring very suggestive 
of the Blazar phenomenon, both qualitatively and quantitatively (amplitude and time scale) f|72jp. These regimes of 
sporadic flaring depend on the power of the turbulence that maintains the in situ Fermi acceleration (of second order 
type). The level of the turbulence energy density varies with the energy mass density of pairs that absorb the power 
input: the increase of the pair energy-mass leads to more absorption of Alfven waves energy and thus reduces the 
efficiency of the in situ acceleration. 
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VIII. DISCUSSION 

In these lectures, the main elements for studying the formation of jets in Black Hole environments have been pre- 
sented together with the important opened questions that deserve continuing the theoretical effort for the achievement 
of this major topic of modern astrophysics. After having exposed the common knowledge in this field that has been 
developed during more than twenty years, some specific views about the way to solve some major difficulties has been 
proposed. Let summary the main problems encountered for relativistic jet formation. 

First, the large scale AGN jets do not seem to be relativistic, especially FR1, and the energy budget of FR2 hot 
spots does not clearly imply a relativistic inflow, whereas their power is comparable to the accretion power, which 
means that they do not lose significant energy during their travel. Unsteady relativistic ejections are obvious at VLBI 
scales, suggesting a flaring activity. Note that the GRBs collimated winds, probably highly relativistic, are also very 
unsteady. 

Second, relativistic MHD flows have not the self-collimating property of non-relativistic ones. An external confine- 
ment is required but cannot be done by the intergalactic medium in the case of the well collimated FR2 jets. 

Third, the power of the Blandford-Znajek mechanism does not seem sufficient to explain the powerful FR2 jets. 
But its power could be sufficient for generating the relativistic clouds at VLBI scales. 

How Black Holes are generated with extreme rotation? Probably not by the accretion process if the Black Hole 
produces jets at the same time ( 34] ). An equilibrium value of the spin parameter should be expected. Black Holes 
merging is probably necessary to get fast rotation. The fact that the Black Hole at the centre of Milky Way is fastly 
rotating is intriguing. Maybe because it was not producing jets... 

Fourth, the generation of a relativistic flow is severely limited by the Compton drag exerted by the photon field 
emitted by the accretion disk. When Compton drag is responsible of its relaxation, within few tens of gravitational 
radii, the jet vanishes and gives rise to hard X-ray emission. This may happen in Seyfert galaxies, where no relativistic 
jets are seen. 

For all these reasons, for AGNs and micro-quasars, we (myself and my collaborators) have proposed to circumvent 
these difficulties in the frame of a two-flow model that combines a non-relativistic or mildly relativistic jet launched 
by the accretion disk with a relativistic flaring ejection that is channelled by the sub-relativistic jet. An observational 
evidence of these combined flows is strongly suggested in Scorpius XI ([2^ H^)- Besides the confinement produced by 
the self-collimating NR-jet, the relations between the two flows are as follows. The NR-jet is more powerful and part 
of its MHD turbulence energises the R-jet plasma. This R-jet plasma is rich of e + e~ pairs; baryons are progressively 
entrained by the relativistic jet and then, after some travel of at least thousand gravitational radii and probably more, 
both jets mix together. As explained in the text, the constant stochastic heating of particles of the relativistic jets 
strengthens the Compton force, both its drag effect and its rocket effect, and then significant Lorentz factors can be 
obtained, whereas, if not heated, the jet rapidly cools and slow down to a modest Lorentz factor. Moreover, this 
stochastic heating makes this ejection highly variable by stimulating the pair creation catastrophe. 

The study of the two-flow model is not yet completed. It would be important to combine the Poynting flux 
propulsion, together with the effect of the Compton force in an anisotropic radiation field, taking account of the 
stochastic heating of pairs by turbulence and of a progressive baryonic pollution. Moreover, we are currently studying 
the changes of regime in the accretion-ejection flows, where a major part of the accretion disk sporadically becomes 
radiatively inefficient (see pT^I. During this stage, because the compactness drops, the pair creation vanishes, making 
the magnetisation higher, and the Compton drag weaker, which allows a relativistic jet generation by Blandford- 
Znajek effect. When, for some reason related to radiative properties of the disk, a SAD is again set up, pair creation 
and Compton drag become stronger again, which bridles the relativistic jet. 

Among the issues that are currently investigated in this field, the magnetic coupling between the black hole and 
the disk (see for instance HJ), and between the black hole, the disk and the jets, are difficult but crucial topics that 
deserve the development of a new generation of MHD numerical codes. 

The generation of relativistic jets in the environment of Black Holes is not only a fascinating issue per se, but is 
also at the origin of high energy phenomena, that are currently at the forefront of the astroparticles physics. 
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